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® Interpolator for and method of interpolating digital samples. 



@ The process of interpolating between two different sampling rates at any point in a sampling interval is 
obtained using a transversal filter arranged as a continuously variable digital delay line in which the tap 
coefficients of the delay line are made to be a function of the coefficients of an nth degree polynomial and the 
delay between the two sampling rates. 
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INTERPOLATOR FOR AND METHOD OF INTERPOLATING DIGITAL SAMPLES 

This invention relates to interpolators for and methods of interpolating digital samples sampled at a first 
rate into digital samples sampled at a second rate. 

In certain digital applications, digital samples of an analog signal which are transferred from a digital 
system operating under a first sampling clock rate to a digital system operating under a second clock rate 

5 requires a so-called interpolator arrangement to interpolate between the sample clock rates. One way of 
providing such interpolation is disclosed in U.S. Patent No. 4,527.020 issued July 2. 1985. to Y.lto. In the Ito 
arrangement, a digital filter, commonly referred to as a transversal filter, is used to perform the interpolation 
in which the sampling interval of a first sampling clock signal (CK2) is divided into a plurality of segments 
each associated with a respective group of filter tap coefficients. Each digital sample to be interpolated at 

10 each occurrence of a second sampling clock signal (CKi) is obtained by supplying to the digital filter the 
group of tap coefficients associated with that segment of the interval which is present at the time of the 
second sampling clock signal. 

The precision of the ITO interpolation arrangement may be improved by dividing the interval of the first 
sampling clock into a large number of segments. However, doing so requires the arrangement to store in 

15 memory a like number of different groups of tap coefficients. The storage of a large number of different tap 
coefficients becomes unwieldy, especially if the digital filter is a multi-tap transversal filter. 

In prior art interpolators used to interpolate between different sampling rate, the interpolated sample is 
generated using one of a plurality of different groups of filler tap coefficients. In accordance with the 
invention, by contrast, we have discovered that a transversal filter arranged as a continuously variable digital 

20 delay line may be used to interpolate between different sampling rates at any point in the timing interval in 
which the coefficients of the filter taps are generated as a function of the coefficients of an nth degree 
polynomial and the delay between the two clock rates. Thus, the coefficients of the delay line required for 
precise interpolation of signal samples at a particular point in the timing interval are generated for that point. 
The invention is. therefore, advantageous since there is no need to store in memory an inordinate number 

25 of different groups of coefficients to achieve precision in the interpolation process. 

The invention will now be described by way of example with reference to the accompanying drawings, 
in which: 

FIG. 1 is a block diagram of a second-order transversal filter embodying the invention; 

FIG. 2 illustrates a signal sampled at a first clock rate In which the samples are interpolated with 
30 respect to a second clock rate: 

FIGS. 3 and 4 illustrate software flow charts depicting the operation of the present invention when 
implemented on, for example, a digital signal processor: and 

FIGS. 5 and 6 each illustrate a matrix of predetermined polynomial coefficients for use in a second- 
order and third-order transversal filter, respectively, embodying the invention. 

3$ 

The response of a finite impulse response filter may be mathematically stated as follows: 

y(t) = I a„x„(t) (1) 
40 n=o 



45 



50 



where N is equal to the number of filter taps. 

Equation (1) is also commonly referred to as being the impulse response of a so-called transversal 
filter, in which the signal elements Xn(t) are available at respective taps of the transversal filter and in which 
the elements an are the values of the respective tap coefficients. 

If the input to the filter is a signal which is a function of time. f(t). and the filter is a delay line, then the 
signal outputted by the filter is a signal which is a function of both time and delay, i.e. f(t - t) where r equal 
to the delay. The output of such a delay line may be generally expressed mathematically as follows: 

n-O 
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With the above in mind. I have recognized that the concepts of transversal filters and delay lines could 
be advantageously extended to a delay line having a variable delay r and tap coefficients which are a 
function of the delay: i.e. Pn{r) in which the values of one such coefficient Pn for various values of r when 
fitted to an nth degree polynomial would be 

5 a function of the delay and respective coefficients of the polynomial. 

Thus, unlike the prior arrangement discussed above which is limited to selecting one of a number of 
different groups of tap coefficients, my arrangement generates the coefficients of the delay line at any point 
in the timing interval based on the value of the delay itself. 

In particular, the tap coefficients Pn (r) shown in equation (2) when fitted to an nth degree polynomial 

iQ may be mathematically stated as follows: pft(T) = an,o + ati.ir ■♦•a„.2'^'»' an.m'''" (3) 

where an.m is a respective coefficient of the polynomial. Applying equation (3) to equation (2) for a matrix of 
polynomial coefficients ranging from ao.o to an.i.m yields the following mathematical expression: 

y(t) = (ao.o + ao.t T + a^^zr^ + ao.mOxo(t) 

+ (a,.o ai.i r + a,.2r2 + a,,„r'")xi (t) 

15 + (a2,o 32.1 r + a2.2r2 + a2.mr'")X2 (t) 

^ (4) 

in which the terms enclosed by the parentheses are the coefficients of the delay line, which are generated 
as a function of the delay. Equation (4) is mathematically equal to the following expression: 

y(t) = (ao.oxo(t) + a^.oxi (t) + a„.oXn(t)) 

20 + (ao.ixo(t) + a,.txt (t) + a„.,x,»(t))T 

+ (ao.2Xo(t) + at.2Xi (t) + a„.2X„(t))r2 

+ (5) 

Equation (5) may be factored to yield a more practical mathematical expression of the output signal as 
follows: 

25 

n-N-1 n-N-l 

a«o D»o D»0 

+ r[ 2 Vm-Xn + ) )) 

n— o 
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Turning now to. FIG. 1. there is shown an illustrative example of a second-order transversal filter 
arranged as a continuously variable digital delay (CVOD) line. The continuously variable digital delay line 
illustrated in FIG. l includes multibit* delay stages, or registers. 201 through 204 for storing a sequence of 
digital samples x(t) through x(t-3T) of the signal received from, for example, a far-end modem. Such a 
signal may be sampled at a rate of 1/T. The latest of such samples, x(t). is stored in register 201, the next 
to the latest. x(t-T). is stored in register 202. and so on. A clock signal, for example, the transmit clock 

^ signal (TO) generated by a modem clock circuit, or the clock signal CK1 disclosed in the aforementioned Ito 
patent, received via lead 213 causes the contents of each register to be shifted from left to right to the next 
register, thereby preparing register 201 for receipt of the next digital sample. 

When digital sample x(t) is received from an external sampling circuit, or input buffer, as the case may 
be, via bus 214, it is stored in register 201 and supplied to multiplier circuits 200-1 through 200-3 via bus 
201-1. Similarly, the digital sample, x(t-T) stored in register 202 is supplied to multiplier circuits 200-4 
through 200-6 via bus 202-1. the digital sample. x(t-2T) stored in register 203 is supplied to multiplier 
circuits 200-7 through 200-9 via bus 203-1 and the digital sample. x(t-3T). stored in register 204 is supplied 
to multiplier circuits 200-10 through 200-12 via bus 204-1. The multipliers 200-1 through 200-12 each 
multiply the digital sample it receives with a respective one of polynomial coefficients ao.o through 33^ and 

50 supplies tile product to a respective summation circuit 205 through 207 via busses 220-1 through 220-12. 
respectively. 

Summation circuits 205 tiirough 207 each sum the digital signals they receive and output digital signals^ 
representative of their respective summations to one of the circuit 209, 21 1 and 212. Multiplier circuit 212 is 
arranged to multiply the summation it receives from circuit 205 via bus 205-1 with the digital value of the 
55 delay received via bus 216 and to supply the product thereof to adder circuit 211 via bus 212-1. Adder 
circuit 21 1 adds the digital value that it receives from multiplier 212 with the summation that it receives from 
summation circuit 206 via bus 206-1 and supplies the sum to multiplier circuit 210 via bus 211-1. Multiplier 
circuit 210 is arranged to multiply the sum it receives from adder 211 and the value of the delay received 
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from delay counter 230 via bus 216. (The manner in which the value of the delay is generated by delay 
counter 230 will be discussed below). Multiplier circuit 210 supplies to adder circuit 209 via bus 210-1 the 
final product involving the value of the delay. Adder circuit 209 adds the value of the product received from 
multiplier 210 with the value of the summation that it receives from summation circuit 207 via bus 207-1 and 
5 outputs a digital signal representative of that sum to bus 215. the digital signal being an interpolated 
sample. 

It is noted that bus 214 could be connected to an input buffer (not shown) which is used to store the 
samples of the received signal. The signal samples may then be clocked one at a time to the input of 
register 201 using the TC clock signal. In addition, the interpolated sample y(t) outputted by adder 209 

to could be stored in an output buffer (not shown) connected to bus 21 5. 

As mentioned above, delay counter 230 shown in FIG. 1 is used to determine the delay, if any. between 
the sampling rates of clock signals TC and RD. Specifically, delay counter 230 comprises a counter for 
counting the rate of the TC clock signal and a counter for counting the rate off the RC clock signal. The 
difference between the contents of each counter is determined each time that the counter counting the rate 

75 of the TC clock signal passes through zero. This difference is then used to generate the value of the delay 
between the rates of the two clock signals. The delay, as determined by counter 230. is then outputted to 
bus 216. 

It is to be understood of course that the value of the delay supplied via bus 216 could also be a fixed 
value contained in a register the output of which is supplied to bus 216. The register could also be under 
20 program control which periodically changes the contents of the register responsive to a particular algorithm. 
Turning now to RG. 2. there is shown an illustrative example of a signal "A" received from, for 
example, a modem, in which the signal "A" is digitally sampled at the rate of the TC clock (i.e., 1/T) to 
provide signal samples, such as. for example, samples x(t) through x(t-3T). When those samples are 
supplied to a CVDO embodying the invention an Interpolated sample y(t) in step with a respective pulse of 
25 the RC clock is generated. It is seen from FIG. 2 that when a dashed line is drawn passing through each of 
the interpolated samples, the resulting signal "B" is virtually identical to signal "A". 

The precision of interpolating the samples of a signal, such as the samples of signal "A" shown in FIG. 
2 using a CVDD filter circuit embodying the invention is commensurate with the degree of the filter. Thus, 
the interpolation process performed by a third-order transversal filter embodying the invention is more 
30 precise than that performed by the second-order transversal filter shown in FIG. 1. Accordingly, in a 
preferred embodiment of the invention, an 8-tap. third-order CVDD transversal filter operating over a 
frequency band from*,0 to 3150 Hz was designed to achieve precision in the interpolation process. Since the 
number of multiplier circuits required to implement such a filter would be unwieldy, the filter was 
implemented on a digital signal processor (DSP), such as. for example, the digital signal processor 
35 designated DSP-20, which is available from AT&T. The DSP-20 is disclosed in THE BELL SYSTEM 
TECHNICAL JOURNAL. September 1981. Vol. 60. No. 7. Part 2. pp. 1431-1462. which is hereby 
incorporated by reference. 

Briefly, in a DSP embodiment of the invention, the eight stages of the delay line are implemented using 
eight memory locations of the DSP's random access memory (RAM). The latest digital sample and the 
40 seven most recent digital samples are stored in the eight memory locations, respectively. A block of RAM is 
also used to store the predetermined magnitudes of the polynomial coefficients for an eight-tap delay line, 
the coefficients being ao.o through a7.3, as will be shown below. A program contained in the DSP's read only 
memory generates the interpolated value using (a) the digital samples stored in the eight RAM locations, (b) 
the values of the coefficients, and (c) the value of the delay. The manner in which the DSP Interpolates 
45 signal samples is similar to that discussed above in connection with FIG. 1. It is noted that the value of the 
delay between clock rates may be determined using the delay counter 230 shown in FIG. i. The generated 
value of the delay Is then supplied to the DSP for use in generating the interpolated signal sample. 

After the interpolated value has been generated and outputted to. for example, an output buffer, a 
program stored in the DSP's ROM effectively shifts the stored digital samples to the right to prepare the 
so first of the eight RAM locations for the storage of the next inputted digital sample. I.e.. the latest digital 
sample. 

We now turn to a discussion of the software programs which implement the mvention on a DSP. 

Turning then to RG. 3, there is shown a flow chart of the DSP program which shifts the inputted signal 
samples stored in respective memory registers to the right in the manner discussed above in connection 
55 with registers 201 through 204 of FIG. 1 to prepare the first of such memory registers for receipt of the next 
digital sample. As noted above, in my DSP implementation, eight taps are used thereby requiring eight 
delay stages, or registers. Thus, the program shown in FIG. 3 causes each stored digital sample to be 
shifted right to the next or succeeding memory register. It is to be understood of course that the program 
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shown in FIG. 3 is not lied to a delay line having eight taps and may be used with a continuously variable 
digital delay line having virtually any number of taps (registers). 

Specifically, when the program illustrated in FIG. 3 Is entered at block 300 it proceeds to block 301 
where it sets the variable n to be equal to one less than the number of taps (M) in the delay line. (It is noted 
5 that the value of M for the delay line illustrated in FIG. i would be 4). The program then proceeds to block 
302 where it begins the process of shifting digital samples from one memory register to the next 
succeeding memory register, beginning with shifting the digital sample from the seventh memory register to 
the eighth memory register. When that shift is completed the program then proceeds to block 303 where it 
decrements the variable n by one and then proceeds to block 304. At block 304. the program tests the 
10 value of n to determine if it has completed the shifting of the digital samples in the manner discussed 
above. The program makes this determination by comparing the value of n with the number one. If n is less 
than one. then the program concludes that it has completed the task of shifting the digital samples and 
proceeds to block 305. Otherwise, the program returns to block 302 to peHorm the next shift. 

At block 305. the program loads the latest digital sample received by the DSP into the first location of 
;5 the sample memory registers and then exits via block 306. 

Turning now to FIG. 4. there is shown a flow chart of the DSP program which generates the Interpolated 
signal sample from the signal samples currently stored in the sample memory registers. In the interest of 
clarity, the following discussion will be directed to a second-order delay line so that the flow of the program 
may follow the flow of the digital samples in the second-order delay line illustrated in FIG. i. 
20 Specifically, the program illustrated in FIG. 4 may be invoked either prior to or after the program 
illustrated in FIG. 3, When invoked at block 400. the program proceeds to block 401 where it clears a 
register designated y which is used to accumulate the results of various operations performed by the 
program. The y register at the completion of the program will then contain the digital value of the 
interpolated signal sample. The program also sets the variable m to equal two to handle the second-order 
25 case and sets the constant T to the number of stages in the filter, which for the circuit shown in FIG 1 is 
four. The program then proceeds to block 402. (It is noted that for the aforementioned 8-tap third-order 
delay line the variable m would be set to three and the constant T would be set to 8). 

At block 402, the program sets the variable n to equal zero and proceeds to block 403. it is noted that 
blocks 403 through 405 comprise a looping routine which performs the function performed by each group of 
30 multipliers and their associated summation circuit based on the values of variables n and m. such as the 
group of multipliers 200-1. 200-4, 200-7 and 200-10 and summation circuit 205 shown in FIG. 1 when the 
values of n and m are 0 and 2. respectively. In particular, during the first pass through the loop, the product 
of the latest signal sample and delay line coefficient ao.2 is added to the contents of the y register, the y 
register being analogous to accumulation circuit 205. The variable n is then incremented by one at block 
35 404 and compared with the constant T at block 405. If the value of n is less than T. i.e. less than four, as 
would be the case after the first pass through the loop, then the program returns to block 403. Othen/vise. it 
proceeds to block 406. During the second pass through block 403. the program adds the product of the 
next-to-the-latest signal sample and coefficient a,. 2 to the contents of register y. This process is repeated 
for coefficients a2.2 and a2.3 and the remaining signal samples. When the last signal sample has been 
40 processed, for example, the sample contained in register 204 of FIG. 1. the variable n is once again 
incremented by one at block 404. thereby making it at least equal to the value of T and causing the 
program to proceed to block 406. 

At block 406 the program multiplies the contents of the y register with the current value of the delay r 
and proceeds to block 407. At block 407. the program decrements the variable m by one and then 
45 proceeds to block 408 where it tests the value of m to determine if it is less than zero. If the value of m is 
not less than zero, then the program transfers to block 402 to perform the function performed by the next 
group of multipliers and their associated summation circuit, i.e. multipliers 200-2, 200-5. 200-8 and 200-1 1 
and summation circuit 206. When the program completes that task, it once again proceeds block 407 via 
block 406 to decrement the variable m and then proceed to block 408. At this point, the program would find 
50 that the value of m is still not less than zero and. therefore, would proceed to block 402. where it would 
perform the function performed by multipliers 200-3. 200-6. 200-9 and 200-9 and summation circuit 207, as 
directed by the values of variables m and n. After performing that function, the program would then proceed 
to block 407 via block 406 and would once again decrement the variable m and then proceed to block 408.. 
At this point, however, the program would find that the value m is less than zero, thereby causing the 
55 program to output the digital value contained in the y register, the digital value being the interpolated signal 
sample. The program then exits via block 409 after completing that task. 

We turn now to discussion of generating the nth degree polynomial coefficients, discussed above, and 
shown in FIG, 1 as ao.o through 33.2- 
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Specifically, the transfer function of a filter with a flat delay r may be stated mathematically as follows: 
G(«.r) - e""**^ (7) 

The transfer function of a transversal filter with coefficients ... C, ... and sampling inten/al T may be 
stated as follows: 

5 

G(a;)- £ C^e'-^^ (8) 

where N is equal to the number of taps of the transversal filter (CVOD). 

Using a delay parameter a such that r = aT and tap coefficients which are polynomials in a then: 

G(..a)-"s"c.(a)..-^-°'2'Eo°C...'""^-«"^''''' W 

where M iis equal to the degree of the polynomial. . 
To determine the coefficients for a CVOD. the following expression (10) derived from equation (9) is 
20 then minimized using the method of least squares with respect to the matrix [C]: 

2j u, a, D m 

Using a CVOO of length N. where N is even, expression (10) is then subjected to the following 
constraints to determine the values of the respective coefficients of the polynomial: 
30 C„(i) = 0:n*^-1 
= l:n= ^ i 
C„(-i)=0: n* N/2 

~ ^ 'Emptoying ixiJession (10) and the constraints imposed by (1 1). the polynomial coefficients ao.o through 
a,, for a second-order four-tap transversal filter (CVOD) having a cut-off frequency of 1350 Hz and an 
interval T equal to 1/9600 Hz were determined as shown in FIG. 5. The polynomial coefficients for a third- 
order 8-tap transversal filter having cut-off frequency of 3000 Hz and an interval T equal to l,«600 Hz were 
also determined and are shown in FIG. 6. ^; 

The foregoing is merely illustrative of the invention, and vanations thereof will occur readily to those 

40 Skilled in the art. 
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1 An interpolator for interpolating digital samples of a signal sampled at a first clock rate into digita 
samples sampled at a second clock rate, including means for storing a plurality of coefficients (aoo to a^a 
each having a predetermined magnitude, for receiving and storing (201-204) a sequence of the digital 
samples (x(t)) sampled at the first clock rate, and for combining.(200-l to 200-12) each one of the sequence 
r^gital samples with respective ones of the coefficients. CHARACTERISED BY means (230) for 
determining the delay between the first and second clock rates, the first-mentioned means serving to 
generate (205 to 207. 209 to 212) the interpolated digital samples as a function of each of the combinations 

and ^^^^^^^^ ^ ^,3.^ j„ g„i„ , therein the coefficients are the coefficients of an nth-degree 

**'Ta method of interpolating digital samples of a signal sampled at a first clock rate into digital samples 
sampled at a second clock rate, including storing a plurality of coefficients each having a predetermined 
magnitude, storing a sequence of the digital samples sampled at the first clock f ^"^i" tSd'bJ 
one of the sequence of the digital samples with respective ones of the coefficients. CHARACTERISED BY 
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determining the delay between the first and second clock rates, and generating the interpolated digital 
samples as a function of each of the combinations and the delay. 
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